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We report the first measurement of the angular coefficients of Z bosons produced in pp col¬ 
lisions and decaying to muon pairs. These coefficients govern the decay of the Z boson and 
thereby the kinematics of the lepton. Their values follow from the vector and axial vector (V-A) 
structure of boson-fermion couplings. The general structure of the lepton angular distribution 
in the boson rest frame is given by 
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d cos 9*d(p* 
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(1 -|- cos^ 6 *) + Ao^{l — 3cos^ 9*) Ai sin(20*) cos^* + Z. 2 ^ sin^ 9* cos(2^*) 
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Here, 9* and (p* are the polar and azimuthal angles of the negatively charged lepton in the rest 
frame of the lepton pair. In this analysis we choose the Collins-Soper (CS) frame [[H to measure 
the angular coefficients A„ considering the momentum of the beam proton closest in rapidity 
to the Z boson as the "target momentum" in [11. The parameters Aq, Ai, and A 2 are related to 
the polarization of the Z boson, whilst A 3 and A 4 are also sensitive to the V-A structure of the 
couplings of the muons. All angular coefficients vanish as the Z boson transverse momentum 
qj approaches zero except for A 4 , which is the electroweak parity violation term. 

The only previous measurement of four of the angular coefficients was performed by the CDF 
Collaboration in pp interactions for qj up to 55GeV |3. The angular coefficients in pp colli¬ 
sions are expected to differ from those in pp collisions for several reasons. For pp collisions, 
the Z boson production occurs predominantly via the qq annihilation process, whilst the con¬ 
tribution of the qg Compton process is larger in pp collisions than pp collisions. Using the 
POWHEG estimation ||3H6l the fraction of qg process in pp collisions at the LHC is 47%; it is 
35% near qj = 0 and increases to ~80% at qj > 100 GeV. For the qq process in the GS frame, 
Aq = A 2 = qj/ (M| -|- qj) llTtiTOl. where Mz is the Z boson mass. For the qg Gompton process 
Aq = A 2 ~ 5qj/ (M| -|- 5qj) IHTH . The relation Aq = A 2 is known as the Lam-Tung relation llT^ . 
reflecting the full transverse polarization of vector boson coupling to quarks, as well as rota¬ 
tional invariance |[l3l . Processes containing non-planar configurations (e.g., from higher order 
multi-gluon emission) smear the transverse polarization, leading to A 2 < Aq fl4l . In contrast 
to what happens at the Tevatron, the average handedness of Z bosons is nonzero at the LHG, 
as for the W boson ||T5l - [T7l . 

The angular coefficients of Z bosons produced in pp collisions at \/s = 8 TeV and decaying to 
are measured as a function of qj and rapidity y. The data, taken with the GMS detector at 
the LHG, corresponds to an integrated luminosity of 19.7fb~^. The large Z boson event sample 
collected by the GMS experiment allows precision measurements of the angular distribution 
for qj < 200GeV and \y\ < 2.1. The coefficients, measured as a function of qj and \y\, are 
compared with three perturbative QGD predictions by FEWZ at next-to-next-to-leading order 
(NNLO) [18], POWEIEG at next-to-leading order (NLO) [SHS], and MadGraph at leading order 

(LO) m- 

The central feature of the GMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and plastic scintillator 
hadron calorimeter, each composed of a barrel and two endcap sections. Muons are measured 
in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. Exten¬ 
sive forward calorimetry complements the coverage provided by the barrel and endcap detec¬ 
tors. Muons are measured in the pseudorapidity range \q\ < 2.4, with detection planes made 
using three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. 
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A more detailed description of the CMS detector, together with a definition of the coordinate 
system and the relevant kinematic variables, can be found in Ref. |20| . 

Matching muons to tracks measured in the silicon tracker results in a relative px resolution for 
muons with 20 < px < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the endcaps. 
A particle-flow (PF) event reconstruction algorithm 11211 l22l is used in this analysis. It consists 
of reconstructing and identifying each single particle with an optimized combination of all 
subdetector information. A trigger for single isolated muon is used, requiring px > 24 GeV 
and Ip I < 2.1. The leading in px reconstructed muon is matched to the muon selected by the 
trigger. 

The signal process is simulated using the MadGraph 1.3.30 generator fT9l with zero to four 
additional jets, interfaced with pythia v6.4.24 |[23l with the Z2* tune 1241 . The matching be¬ 
tween the matrix element calculation and the parton shower is performed with the kx-MLM 
algorithm Il25l . The GTEQ6L1 Ii26l parton distribution functions (PDF) are used for the event 
generation. Multiple-parton interactions are simulated by pythia. The powheg generator |3t- 
(6l interfaced with PYTHIA (same version used for MadGraph) and the GTIO PDF set |27l are 
used as an alternate to test any model dependence in the shapes of the angular distributions. 

Background simulations are performed with MadGraph (W+jets, tf, tt), POWHEG (single 
top quark Il28l l29l l. and PYTHIA (WW, WZ, ZZ). The normalizations of the inclusive Drell- 
Yan, W boson [TSl, and tt II30I distributions are set using NNFO cross sections. For single top 
quark production a higher order (approximate NNFO liSTl ) inclusive cross section is used. The 
generated events are passed through a detector simulation based on Geant4 1321 . 

Fach muon candidate is required to be reconstructed in the muon detectors and in the inner 
tracker, and the global track fit is required to have a reduced < 10- The vertex with the 
highest sum of pj for associated tracks is defined as the primary vertex. The distance of the 
muon candidate trajectories with respect to the primary vertex must be smaller than 2 mm in 
the transverse plane and 5 mm along the beam axis. The leading (subleading) muon is required 
to have pj > 25 (10) GeV and |;^| < 2.1 (2.4). In order to suppress background events, the 
muons are required to be isolated from nearby particles. The relative isolation is calculated as 
the ratio of the scalar sum of px of all PF candidates from the same primary vertex, within a 
cone of AR = \/ (Ap)^ -|- {Acp)^ < 0.4, and the pj of the muon. For the leading (subleading) 
muon in px, the relative isolation must be less than 0.12 (0.5). Oppositely charged muon pairs 
with an invariant mass in the range 81-101 GeV are selected. In the rare case that more than two 
muons are selected, the muon pair with invariant mass closest to the Z boson mass is chosen. 
The muon pair must satisfy \y\ < 2.1 since at higher \y\ the acceptance varies rapidly. After 
the event selection, 4.3x10^ events with Z boson candidates remain for \y\ < 1.0 and 2.5x10^ 
events for 1.0 < \y\ < 2.1. 

A "tag-and-probe" method ||33| is used to measure the efficiencies for track reconstruction, 
trigger, muon isolation, and muon identification in data and simulation. Efficiency corrections 
are applied as multiplicative scale factors to the simulation values. The efficiency for track 
reconstruction is measured in bins of p since the px dependence is weak. The trigger efficiency 
is determined in bins of px and p, separately for p+ and p^. The identification efficiency is 
measured in bins of px and p. Since the subleading muon can point in the direction of the 
hadronic activity, a looser isolation requirement is used and its efficiency is measured as a 
function of px/ cos 0*, and (p*. The efficiency of the isolation requirement for the leading muon 
is measured as a function of px and p of the muon, as detector effects relate to these variables 
more directly than to the Z boson cjj and y. 
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After event selection, the background contribution ranges from ~0.1% at low qj to ~1.5% at 
high qj. The yields of the backgrounds from tf, tt, WW, tW, and W+jets production are esti¬ 
mated from data using lepton flavor universality. Most of these backgrounds typically have 
two prompt leptons, which may have the same flavor. The W+jets background is flavor asym¬ 
metric, but its contribution is small. We assume that the ratio of the number of oppositely 
charged background fiji and e}i events is the same in data and simulation. We use the ratio of 
the e}i yields in data and simulation after applying muon and electron selection criteria ll33llMI 
to normalize the simulation to data. 

The acceptance and the efficiency at the event level vary in cos 6 * and (p*, and strongly with qj 
and y. In order to avoid a bias in the acceptance due to the modeling of the Z boson kinematics, 
the simulation is reweighted in fine bins of qj and y to match the background-subtracted data 
distribution. The weights are determined at the reconstruction level and applied at the genera¬ 
tor level. The weighting is iterated four times, with negligible change between the second and 
fourth iteration. 

The angular coefficients are measured in eight bins of qj and two bins of \y\, by fitting the 
two-dimensional (cos 6 *, (p*) distribution in data with a linear combination of templates. These 
templates are built for each coefficient A, by reweighting the simulation at generator level to 
the corresponding angular distribution, as given in Eq. Q. The templates are based on recon¬ 
structed muons, and thereby incorporate the effects of resolution, efficiency and acceptance. 
A template is also built for the term (1 -|- cos^0*) of Eq. 0 . An additional template, with 
shape and normalization fixed, is developed for fitting the backgrounds. A binned maximum- 
likelihood method with Poisson uncertainties is employed for the fit. The angular coefficients 
As, Ag, and A 7 are predicted to be very small; they are set to zero and excluded from the fit. 
Since Aq through A 4 are sign invariant in cp*, the absolute value |^* | is used. The fit is made in 
12 X12 equidistant bins in cos 6 * and \(p*\- The statistical uncertainties from the fit are confirmed 
by comparison with pseudo-experiments. 

To test the robustness of the result with respect to the analysis method and trigger effect, the 
angular coefficients Aq, A 2 , A 3 , and A 4 are also measured by an independent analysis similar 
to that reported in Ref. f2l, where one-dimensional (ID) templates produced using powheg 
are fitted to the distributions in cos 9* and |(/>* |. The ID fit analysis is performed iteratively, so 
as to be unbiased with respect to the assumed templates and to possible correlations between 
COS0* and |^*|. The analysis differs in the triggers, estimation of backgrounds, simulation, 
and selection criteria. The ID fit analysis uses a sample that requires a dimuon trigger with 
asymmetric muon pj thresholds of 17 and 8 GeV. Both results are consistent within their total 
systematic uncertainties, excluding uncertainties common to both analyses. 

Some examples of the measured cos 6 * and |^* | distributions from the ID analysis are given in 
Kg|Il The measured and simulated distributions are shown together using the best fit values 
of the angular coefficients. The shape of the cos 6 * distribution changes with qj and \y \ because 
the acceptance and efficiency in cos0* depend strongly on these two variables. Eor |^*|, the 
shape of the distribution changes moderately with qi, and is almost insensitive to \y\. The 
comparison of data and simulation shown in Eig. gives confidence that the acceptance and 
efficiency dependences are correctly modeled in the simulation. 

Several sources of systematic uncertainties are taken into account. The most significant source 
is the muon efficiency that includes the trigger, track reconstruction, isolation, and identifica¬ 
tion. The statistical uncertainties of the measured efficiency scale factors are taken into account 
by simulating 500 pseudo-experiments in which the templates are reformed, each time vary¬ 
ing the scale factors randomly within the given uncertainty. The systematic uncertainties in 
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Figure 1: A few examples of the observed ID angular distributions in cos0* (left) and |^*| 
(right) compared to the MC simulation using the best fit values of the angular coefficients. The 
top (bottom) plots show the distributions for 10 < qj < 20 GeV (120 < qj < 200 GeV), a region 
where Aq and A 2 are small (large). The background-subtracted data points are shown with 
filled (open) circles for \y\ <1 {1 <\y\ < 2.1), whilst the corresponding MG results are shown 
with the solid (dashed) lines. Vertical bars represent the statistical uncertainties. The lower 
panels show the data-to-MG ratios. 


the extraction of the efficiency (e.g., background estimates) are also included. Another sig¬ 
nificant uncertainty stems from the statistical precision of the templates, which is estimated 
using pseudo-experiments. The pileup uncertainty is estimated by varying the cross section 
of the minimum bias events by ±5%. The muon momentum bias is measured in data and 
simulation, and corresponding corrections are applied Il35ll . The statistical uncertainties in the 
muon momentum correction factors are propagated to a systematic uncertainty using pseudo¬ 
experiments. In addition, a systematic uncertainty is assessed to take into account possible 
global offsets from the peak position of the Z boson mass. The systematic uncertainties for the 
background are estimated by varying the normalization scale factor of the ey sample by 10% 
and the yields of WZ and ZZ events by 50%. The statistical precision of the iterative reweight¬ 
ing is determined using pseudo-experiments. The difference between the last two iterations 
is assigned as additional systematic uncertainty. The effect of final-state radiation is taken 
into account by adding the energy of photons within a cone of radius 0.1 around the muon 
direction ||3^ . Weights are applied to the simulation to reflect the difference between a soft- 
collinear approach and the exact O(aQED) result and the reconstructed template is rebuilt using 
the weighted simulation. The difference between templates is used to estimate the systematic 
uncertainty from final-state radiation. Finally, the acceptance uncertainty, related to the values 
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of Ai assumed in the simulation, is estimated by reweighting with the fitted values of and 
the difference in results is included as a systematic uncertainty. Generally, the statistical uncer¬ 
tainties dominate in the highest bins in qj, whilst the systematic uncertainty in the efficiency 
tends to be the most important elsewhere. 

The results of the qj and \y\ dependent measurements of the angular coefficients Aq to A4 as 
well as the difference Aq — A^ are presented along with MadGraph, powheg, and FEWZ (at 
NNLO) calculations in Figs. and The various systematic uncertainties of the five angular 
coefficients Aq to A4 are presented in Fig.|^ The values and uncertainties of the coefficients are 
provided in Tables [^and|^ The PDF sets used in the calculations are GTEQ 6 L for MadGraph 
and GTIO for POWHEG (at NLO) and FEWZ (at NNLO). The MadGraph predictions for A4 
are systematically higher than those of POWHEG and FEWZ because MadGraph uses a weak 
mixing angle calculated without considering radiative corrections. The measured Aq and A 2 
coefficients agree better with the prediction of MadGraph than with those of POWHEG and 
FEWZ, especially at high qj. At qj = 0, the POWHEG prediction for Aq is negative, which is 
unphysical and has been traced to approximations in the shower matching algorithm. The 
FEWZ prediction is shown for qj > 20 GeV, where the calculations are considered reliable. We 
find that Aoiqj) and A 2 iqT) are larger in pp collisions than those measured in pp collisions at the 
Tevatron. The larger contribution from the qg process in pp collisions at the LHG is responsible 
for this difference. We observe the violation of the Lam-Tung relation {Aq = A 2 ) anticipated by 
QGD calculations beyond leading order Il37l . We find that Aq > A 2 , especially for high qj. In 
addition, we measure nonzero values of Ai and A 3 . The comparison of the results for \y\ < \ 
and 1 <\y\ < 2.1 is shown in Fig.|^ 

In summary, we presented the five major angular coefficients, Aq through A 4 , for the produc¬ 
tion of the Z boson decaying to muon pairs as a function of qj and \y\ in pp collisions. These 
results play an important role in future high-precision measurements, such as the measurement 
of the mass of the W boson and of the electroweak mixing angle. Some theoretical predictions 
deviate from the measurements in qj. Further refinements of the theory are needed to achieve 
a better agreement with the experimental results. 
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Table 1: The five angular coefficients Aq to A 4 and Aq — A 2 in bins of qj for \y\ < 1. 


qj [GeV] 

^0 

ii^stat 

i^syst 

Ai 

±‘^stat 

i^syst 

A 2 

±<^stat 

i^^syst 

0 - 

TO 

0.018 

±0.003 

±0.009 

-0.008 

± 0.002 

±0.005 

0.007 

±0.004 

±0.003 

10 - 

-20 

0.068 

±0.004 

± 0.010 

-0.006 

±0.003 

±0.005 

0.037 

±0.004 

±0.005 

20 - 

-35 

0.179 

±0.004 

±0.013 

0.014 

±0.003 

±0.008 

0.136 

±0.006 

±0.014 

35- 

-55 

0.357 

±0.006 

±0.013 

0.033 

±0.005 

±0.014 

0.278 

±0.008 

± 0.022 

55- 

-80 

0.563 

±0.007 

± 0.010 

0.031 

±0.007 

±0.017 

0.447 

± 0.012 

± 0.022 

00 

0 

1 

■120 

0.716 

± 0.010 

±0.009 

0.029 

± 0.010 

±0.017 

0.583 

±0.017 

±0.037 

120 - 

-200 

0.834 

±0.015 

±0.014 

0.002 

±0.015 

±0.013 

0.741 

±0.029 

±0.043 

>200 

0.928 

±0.035 

±0.015 

- 0.020 

±0.032 

± 0.012 

0.689 

±0.068 

±0.035 

qj [GeV] 

^3 

±‘^stat 

i^syst 

A 4 

±‘^stat 

i^syst 

Aq — A 2 

±^stat 

±^syst 

0 - 

TO 

0.007 

± 0.002 

±0.004 

0.020 

± 0.002 

± 0.002 

0.011 

±0.005 

±0.009 

10 - 

-20 

0.003 

± 0.002 

±0.003 

0.013 

±0.003 

± 0.002 

0.032 

±0.006 

± 0.011 

20 - 

-35 

0.006 

±0.003 

±0.003 

0.015 

±0.003 

±0.003 

0.043 

±0.007 

±0.016 

35- 

-55 

0.005 

±0.004 

±0.005 

0.021 

±0.004 

±0.004 

0.079 

± 0.010 

±0.018 

55- 

-80 

0.009 

±0.006 

±0.006 

0.002 

±0.006 

±0.004 

0.116 

±0.014 

± 0.022 

00 

0 

1 

■120 

0.033 

±0.008 

± 0.010 

0.019 

±0.008 

±0.005 

0.133 

±0.019 

±0.032 

120 - 

-200 

0.008 

±0.014 

± 0.010 

0.010 

± 0.012 

±0.007 

0.093 

±0.031 

±0.031 

>200 

0.101 

±0.031 

±0.016 

0.029 

±0.026 

± 0.010 

0.239 

±0.072 

±0.028 


Table 2: The five angular coefficients Aq to A 4 and Aq — A 2 in bins of qj for 1 < \y\ < 2 . 1 . 


qj [GeV] 

^0 

±‘^stat 

i^syst 

Ai 

±‘^stat 

i^syst 

A 2 

±^stat 

i^syst 

0 - 

TO 

0.032 

±0.005 

± 0.010 

0.002 

±0.003 

±0.007 

0.019 

±0.005 

±0.006 

10 - 

-20 

0.077 

±0.006 

±0.009 

0.018 

±0.004 

±0.006 

0.038 

±0.005 

±0.007 

20 - 

-35 

0.179 

±0.008 

±0.013 

0.038 

±0.005 

±0.008 

0.129 

±0.006 

±0.016 

35- 

-55 

0.385 

± 0.011 

±0.017 

0.063 

±0.007 

± 0.011 

0.260 

±0.009 

±0.024 

55- 

-80 

0.554 

±0.013 

±0.015 

0.066 

± 0.011 

±0.016 

0.448 

±0.014 

± 0.021 

00 

0 

1 

■120 

0.737 

±0.015 

±0.014 

0.059 

±0.015 

±0.019 

0.587 

± 0.021 

±0.031 

120 - 

-200 

0.860 

± 0.020 

± 0.012 

0.064 

± 0.021 

±0.018 

0.758 

±0.035 

±0.035 

>200 

0.876 

±0.045 

± 0.020 

0.040 

±0.044 

± 0.020 

0.864 

±0.087 

±0.041 

qj [GeV] 

^3 

±‘istat 

±^syst 

A 4 

±^istat 

±^syst 

Aq — A 2 

±<^stat 

±^syst 

0 - 

TO 

0.009 

± 0.002 

±0.005 

0.076 

±0.003 

±0.004 

0.013 

±0.007 

± 0.011 

10 - 

-20 

0.003 

± 0.002 

±0.004 

0.072 

±0.004 

±0.005 

0.039 

±0.008 

± 0.011 

20 - 

-35 

0.012 

±0.003 

±0.006 

0.044 

±0.005 

±0.007 

0.051 

± 0.010 

±0.017 

35- 

-55 

0.012 

±0.005 

±0.008 

0.052 

±0.007 

±0.009 

0.124 

±0.014 

± 0.021 

55- 

-80 

0.036 

±0.007 

±0.018 

0.052 

±0.009 

±0.008 

0.106 

±0.019 

±0.019 

00 

0 

1 

■120 

0.074 

± 0.010 

±0.028 

0.074 

± 0.011 

±0.014 

0.150 

±0.025 

±0.028 

120 - 

-200 

0.121 

±0.017 

±0.029 

0.056 

±0.016 

±0.017 

0.102 

±0.039 

±0.031 

>200 

0.181 

±0.041 

±0.027 

0.005 

±0.034 

±0.017 

0.012 

±0.090 

±0.039 
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Figure 2: Comparison of the five angular coefficients A; and Aq — A 2 measured in the Collins- 
Soper frame in bins of qj for \y\ < 1. The circles show the measured results. The vertical bars 
represent the statistical uncertainties and the boxes the systematic uncertainties of the measure¬ 
ment. The triangles show the predictions from MadGraph, the diamonds from POWHEG, and 
the crosses from FEWZ at NNLO. The boxes at the FEWZ values indicate the PDF uncertainties. 
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Figure 3: Comparison of the five angular coefficients and Aq — A 2 under the same conditions 
as Fig. 1^ for the rapidity bin 1 < \y\ <2.1. 
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